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1 Introduction 



The averaging problem in general relativity was first raised by Shirkov and Fisher [T] 
in 1963. The energy-momentum-tensor used in cosmological models is an average over 
the non homogeneous tensor present in nature. But, due to the non linear nature of 
Einstein's equations, the metric belonging to the averaged energy-momentum tensor is 
not identical to the averaged metric. The averaging prescription by itself provides a 
fundamental problem because of the freedom of choice of coordinates. The authors of [I] 
suggested to integrate the metric tensor over a four dimensional volume with the familiar 
factor sj—g in the measure. Such an expression is, however, not covariant due to the 
freedom of performing local transformations. A covariant averaging prescription can 
be constructed by introducing a bivector gP(x,x') of geodesic parallel displacement, as 
discussed in the appendix of [2]. This transforms as a vector with respect to coordinate 
transformations at either x or x' and maps a vector Ap{xf) to A a (x) = gP(x,x')Ap(xf), 
analogously for higher order tensors. An averaging with the help of bivectors was also 
used in the work of Zalaletdinov [3] where the emphasis was on the commutativity of 
averaging and covariant differentiation. As remarked by Stoeger, Helmi, and Torres [I] 
the method of using a covariantly conserved bivector is not applicable to the metric, 
because the covariant derivative of the metric vanishes. The metric is therefore invariant 
under this averaging procedure. In the thesis of Behrend [S] the metric is represented by 
tetrads and the averaging performed over the latter. The tedrads are chosen according 
to a covariant minimalization prescription. 

This is only an extremely brief survey of the literature. For more references as well 
as the implications for the fitting problem, back reaction, contributions to dark energy, 
we refer e.g. to the monograph of Krasinski [B] and to the comprehensive recent review 
of Buchert [7] . 

Under a covariant averaging process we understand a prescription which has the 
following properties. Let two observers describe the same physics in different coordinate 
systems S and 5", with metric tensors g^ u and g'. Both of them apply a definite 
averaging procedure in their respective systems, resulting in the averaged metrics 
< g^v > an d < g\ iv >, respectively. Then the results have to be connected by the same 
transformation as the original metric, i. e. 



In other words, the operations of averaging and of coordinate transformations have to 
commute. Furthermore, averaging over a region which is closely located around some 
point should, of course, reproduce the metric at this point. 

In two respects the approach in the present paper is modest. Firstly we assume 
that a reasonable foliation into space and time has already been performed, so we will 
essentially concentrate on spatial averaging of a perturbed flat metric, with time kept 
fixed. At the end we give a simple special generalization to four dimensions. Secondly 
we will demonstrate covariance only in first order of the perturbation of the metric. 

In other respects our approach is ambitious. We give a closed formula for the av- 
eraging prescription, there is no need to solve any differential equations or to resort 
to a background of dust, perfect fluid, or whatsoever. We not only can average over 
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a given sphere, but may introduce an arbitrary smearing function f(r). This appears 
more physical than a sharp cutoff at the boundary of the considered sphere. A very 
important result, at present only derived in three dimensions, is the fact that we can 
apply the same formula which works for the metric as well to the Einstein tensor. This 
implies that the Einstein equations for the averaged metric are identical to the averaged 
equations. 

The mapping of the metric gki(x r ) to the averaged metric < g mn > (x) is represented 
by a bitensor K® n (x' — x) which we will specify in detail. A product of bivectors, as 
frequently used in the literature, is not sufficient. A mapping with the help of a bitensor 
was also formulated by Boersma [8], although without going into details. 

The paper is organized as follows. 

In sect. 2 we introduce the necessary technicalities for the projective coordinates. 
Technically these are much more convenient than the familiar polar coordinates, in par- 
ticular for the many partial integrations which we have to perform. In sect. 3 we present 
the general form of the averaging formula. A central aspect is that, besides the arbitrary 
smearing function f(r), the integrals F(r) and G(r) over / and f/r appear. This is 
essential for the proof of covariance in sect. 4. Fixing the remaining freedom of the 
parameters appropriately, we show in sect. 5 the surprising and highly welcome result, 
that our formula also yields a covariant averaging of the Einstein tensor. In sect. 6 we 
discuss a simple four dimensional covariant generalization for static perturbations. The 
proof for the covariant averaging of the Einstein tensor "almost" goes through also in 
this case, but fails at the very end. Principal limitations and possible generalizations are 
discussed in the conclusions. 



2 Projective coordinates 

For the moment we choose the point, where the averaging is to be performed, as the 
origin for simplicity. Instead of using standard polar coordinates x = r sin 9 cos (ft, y = 
rsin#sin0, z = rcos8, it is much more convenient for various reasons to use the coor- 
dinates r,u = v} , v = u 2 of stereographic projection. They are connected to the polar 
coordinates by u = tan | cos (ft, v — tan | sin (ft, and to the cartesian coordinates by 



30 y I 

u = u 1 = , v = u 2 = , r = Jx 2 + y 2 + z 2 , (2.T 

r + z r + z v 



or, vice versa, 



2U 2V 1 — W . . 2 2 2 /n 

x = r, y = r, z = r, with w = u + v . (2.21 

1 + w 2 1 + w l 1 + w l 



The geometric meaning is simple. 2yw 2 is the distance from the north pole (0,0, 1) to 
the point where the straight line from the south pole (0, 0, —1) through (x, y, z) on the 
unit sphere cuts the plane z — 1. The north pole corresponds to u — v — 0, the south 
pole to u 2 -\- v 2 = oo. The integration element for the angular averaging is 

dQ, 1 „ ,„ , , 1 dudv 

— = — S m6d9d(ft = -- ^. 2.3 

47r Air ir{l+w 2 } 2 
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The partial derivatives in projective coordinates become 



d x\ d du l d 
dx l r dr dx l du l ' 



with 



M = «»^-ypfr + ,)^ (for i= 1,2),M = (2.5) 
ox 1 ryr + zy ox 6 r 

Fortunately (2.5) will not be needed explicitely, we will only need two simple obvious 
properties of the partial derivatives du % jdx 1 : They are orthogonal to x l , i.e. x l du l jdx 1 = 
0, and proportional to 1/r for fixed u, v. 

The perturbed three dimensional metric in cartesian coordinates is written as g mn (x) = 
brnn + h mn (x). This decomposition is preserved under translations, rotations, and in- 
finitesimal transformations. In the projective coordinates we also split off the flat part 
and write (indices i,j run from 1 to 2 in (2.6) and refer to u 1 ,^) 

2r Ar 2 
9rr = 1 4~ h rr , g r i = - ■ - h r i, Qij — ■ rrr- + hij\. (2-6) 
1 + -ur [1 + w z \ z 

A corresponding decomposition can be written down for the polar coordinates. 



3 General form of the averaging formula 

The formula presented below looks somewhat strange at first sight, therefore it is appro- 
priate to motivate it. Originally we tried to transform to a specific system, perform the 
average there, and transform back to the original one. In this way covariance is achieved. 
The conditions which essentially fixed the specific system, denoted by primes, were 

Kr - KM = h' rl - h' rl (o) = h' r2 - h' r2 (o) = o. (3.i) 

To avoid singularities of the transformation between the systems, it was necessary to 
split off the values at the origin in these conditions. One could next determine the 
transformation leading from the original to the primed system, perform the average there, 
and transform back. The fact that we had to split off the values at the origin had the 
unwanted consequence that part of the metric remained unaveraged. Instead of pursuing 
this approach further it is, however, more useful to abstract from this original motivation 
and to concentrate on the structure which was obtained in this way. The important point 
is the following. After determining the transformation, choosing a smearing function 
f(r), and performing the average, not only the original function /(r), but also two 
integrals F(r) and G(r) appear. This will be the key for achieving covariance. So much 
for the motivation, all the following is independent of it. 
In detail, we choose a smearing function /(r) with 

[°°f(r)dr = l. (3.2) 
Jo 

Because we are in three space dimensions one should have /(r) ~ r 2 for r — > 0. 
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We will further need the function F(r), the integral over f(r), as well as the function 
G(r), the integral over f(r)/r, both normalized such that they vanish at infinity: 

F(r) = - / f(r')dr', G(r) = - / l^-dr'. (3.3) 

Jr Jr T 

Obviously one has F(0) — — 1. 

The general form of the three dimensional averaging formula now reads 



< 9mn > (0) = J K% n g kl (x, y, z)dr^, (3.4) 



dtt 

< 9mn > (UJ = J K™ n g kl (x, y, z)dr 

with 



Kmn — K[f]mnf( r ) + ^[F/r]mn ~ + ^[G]mn^( r )' (3-5) 

The structure of the tensors (tensors in the sense of linear algebra) K^ mn , (a = /, Fj r, G) 
is, with constant coefficients A^, • • • , F[ ], 



/ h 1 h hi 

+£ [a] [5 m — + 5 m — + <5 n — + tn—^-} + F [a] - . (3.6) 

This is the most general tensor structure which is symmetric under the exchange m <-> n, 
and under k <-> Z. Naive averaging would correspond to = 1/2 and all the other 17 
coefficients , • • • , vanishing. 

There are some restrictions for the coefficients from the beginning. The function 
F(r)/r is singular at r = 0. The expansion of g k \ around r = starts with the constant 
gki(0). For this term we can perform the angular averaging using 

ry ry* f ry* ry* ry ^ ry ^ j 

^-^n», ^^^^(S m nS kl + 5 k m 8 l n + 5 k JJ etc. (3.7) 

The tensor Krp, r -, mn has to vanish after angular integration in order to avoid a singularity 
at r = 0. Collecting the factors in front of 5 k rn 5 l n + S k 5 l m , and those in front of S mn 8 kl , this 
leads to the two conditions 



2 1 

A [F /r] + g#[F/r] + T^I^Al = °> ( 3 ' 8 ) 

+ gC[F/V] + g^IF/r] + — F [ F M = °' ( 3 ' 9 ) 

There are also restrictions for the coefficients which multiply f(r) and G(r). These will 
become relevant for proving covariance. 

K[f]mn * s restricted by the transversality condition xiK^ mn = , which implies 

C lf] = ~ B [fh E lf] = ~ A [fh F lf] = 2A lf] ~ D lf}- ( 3 - 10 ) 
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Finally, K k G ^ rnn contains only terms ~ x k x l , which means that 



k^l 



Aye] = B\a\ = D\a\ = E\ G = 0. 



>[G] 



[G] 



J [G] 



Therefore the expressions for K k ^ mn and K k G]mn simplify to 



' [G] mn 



jy-kl 
A [/]mn 



A 



If] 



(5 k - 



Xm.X 



x n x 



+(B[f]S mn + D [f] 



kl 



+ (m <r 
2 )> 



n) 



K 



ki 



x k x l 



[G\mn 



= K{G]mn With K[ G ] mn = C[ G ](5 m „ + F[ G ] 



X>rnXr 



(3.11) 



(3.12) 



(3.13) 



Let us now investigate the limit of a smearing function f(r) which is closely local- 
ized around r = 0. The same then holds for F(r) and G(r). We may thus put 
gti{x, y, z) = gki{0) in (3.4) and take it out in front of the integral. The angular averages 
can be performed using again (3.7). One is left with the radial integrals J °° f(r)dr = 1, 
and J °° 1 • G(r)dr = — J °° r ■ G'{r)dr = —1. The term with F(r)/r does not enter, be- 
cause the angular integration vanishes due to the conditions (3.8), (3.9) derived before. 
The result has to be identical to g mn (0). This leads to two further conditions, derived 
from comparing the terms with g m „(0) and those with 5 mn <?j(0). Making use of the sim- 
plifications which arise from the restrictions for the coefficients refering to / and G one 
obtains 



-[UA [f] - 2D [f] - 2F [G] ] = 1, (3.14) 
l[2A [f] + 10B [f] +AD [f] -5C [G] -F [G] ] = 0. (3.15) 

The considerations above also show that the averaging of a constant (e.g. of 8 mn ) gives 
back this constant. This also implies, that (3.4) can as well be applied to the perturba- 
tion, i.e. one can replace g mn -> h mn , g M -> h M there. 



4 Conditions for covariance 



Let us apply an arbitrary infinitesimal transformation x k = x' k + £ k , which leads to a 
change Sgki = £k,i + ~~ * 2^,/, when contracted with the symmetrical tensors in (3.6). 
Because invariance with respect to translations and to rigid rotations around the point 
of consideration is manifest, one can restrict to transformations which leave the origin 
fixed, i.e. £fc(0) = 0, such that 

Ux)=ZkMx l + 0(r 2 ). (4.1) 

Consider the change of the integrand in (3.4). In a first step we transform the terms 
~ f(r) and ~ G(r) by partial integration with respect to r, such that all three contri- 
butions become ~ F(r)/r. Because all the tensors K k ^ mn in (3.6) depend on u, v only, 
but are independent of r, this is rather simple. 
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r°° c°° Fir) d 

I /(r)&A = &,i(0) - Jf -Li(r^w)dr, (4.2) 

where we used -F(O) = —1 in the boundary term. 

In the term with G(r) we only use the factor x l /r in (3.13), not the x k jr. Both 
of them are well defined and independent of r, they only depend on the angles (or on 
u, v, respectively) which are still fixed here. One can apply the following chain of partial 
integrations, where we used (4.1) and F(0) = —1 in the boundary terms: 

/ G{r)-i Kl dr = / G{r){-i k )dr = - ^ k dr 
Jo r Jo or Jo r 

= --&,i(0)+r^(|-&-^)dr. (4.3) 
r Jo r or r 

All together this leads to the integrand (to be averaged over the angles) 

2 ( K [f]mn ~ K [G \ mn —^-)t,k,l{ti) 

+ 2 J ^{[- K if]rnnr ¥r +K^ r]m ^ (4.4) 

All manipulations which involve partial integrations with respect to r have now been per- 
formed. We next consider the curly bracket in (4.4) which has to vanish after averaging 
over u, v. We introduce the projective coordinates and use 



-V • 777— 6- (4.5) 



r dr dx l du 

From (2.5) we recall that the partial derivatives du l /dx l are orthogonal to x l and pro- 
portional to 1/r. Moving rd/dr to the right of du 1 /dx l in the term —K^ mn (rd/dr)^k,i 
thus gives an extra contribution. The curly bracket in (4.4) becomes 



U-K* ^ + ^ir-+^ — — 

\l IflmnK r g r -T Qxl dui ) ^ T [f]mn dxl 

The first term contains a second derivative with respect to r. This term has no chance 
to cancel against anything else, but it vanishes due to the transversality condition 
xiK^i mn = 0. We are left with terms ~ d/dr and those ~ 1/r (recall that du % /dx l ~ 1/r). 
Both of them have to vanish after angular averaging for any Thus we multiply by 
the integration element (2.3) and remove all partial derivatives d/du l acting on £ k by 
partial integrations. This gives two conditions which arise from collecting terms ~ d/dr 
and terms ~ 1/r: 
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d f K[f]mn 9u l \ K [F/r]mn Xi K [G]mn X k 

duAll + w 2 }^ dx l ) [1 + w 2 ] 2 r [1 + w 2 ] 2 r 

9 f K [f]mn + K [F/r]mn 9u l \ _ K [G]mn X k 

dtfx [1 + -U7 2 ] 2 dx l J [1 + w 2 ] 2 r 
The sum gives a condition for if^ r ] ran alone, 

9 / K [F/r]mn r 9u l \ K [F/r]mn Xj _ Q 

du % V [1 + w 2 ] 2 dx l J [1 + w 2 ] 2 r 

therefore we will discuss (4.9) and (4.7) in the following. 

One now has to use the properties of the partial derivatives du % /dx l in (2.5), and 
to introduce the expressions (3.6), (3.12), (3.13) for the K^ ]mn . The further treatment 
can be greatly simplified by making use of the fact that, from rotation invariance and 
symmetry, the terms in (4.7) - (4.9) must be a superposition of the form a(r)S mn x k + 
b(r)(5^ n x n + 5 k x m )+c(r)x m x n x k /r 2 . Therefore there are only three invariants which have 
to vanish, and we are free to choose simple special cases in order to determine them. 

The first choice is to take the trace m = n and put k — 3. This gives [3a(r) + 
2b(r) + c(r-)]z. Next one can put m = n = k = 3 which gives [a(r) + 2b(r) + c(r)z 2 /r 2 ]z. 
The vanishing of these two expressions implies already the three equations a(r) = b(r) = 
c(r) = 0. While the first choice gives only one relation, the second one gives two relations 
from collecting the terms ~ (1 — w 2 )/[l + w 2 } 3 and ~ (1 — w 2 ) 3 /[l + w 2 } 5 . The three 
relations such obtained from (4.9) read 

6A[ F / r ] + 9B[ F / r ] + 3C[ F /r] + 3-D[F/r] + 4-E^/V] + F[F/r] = 0, (4-10) 
QA[ F /r] + 3B[F/ r ] + C[F/ r ] ~ 2D[ F /r] ~ 2E[ F / r ] = 0, (4-11) 

5D [F/r] + 10E [F/r] + F [F/r] = 0. (4.12) 

The earlier equations (3.8), (3.9) are consequences of (4.10) - (4.12), therefore we can 
forget them. The same procedure can be applied to (4.7) and leads to three further 
independent conditions: 

4:A[f] + QB[f] + 2D[f] — 2A[ F /r] — 3B[ F / r ] — 3C[f/V] — D[ F / r \ 

—AE[ F / r ] — F[ F / r ] — 3C[ G ] — F[G] = 0, (4.13) 

8 A/1 + 2-B[/] — 2D[f] — 2A[F/ r ] — B[ F /r] — C[F/ r ] ~ 2E[ F /r\ ~ C[G] = 0, (4.14) 

SA [f] -4:D [f] + D [F/r] + 2E [F/r] + F [F/r] + F [G] = 0. (4.15) 

If (4.10) - (4.15) are fulfilled, the integral in (4.4), when averaged over the angles, re- 
spectively over u, v, vanishes, i.e. the expression is gauge invariant. 

Finally we have to consider the boundary term 2(KM, mn — K[ G ] mn x k x l /r 2 )£fc 5 z(0) in 
(4.4). Because £fc,z(0) is constant, the angular averaging can be performed explicitly. In 
order to fulfill the covariance condition (1.1), the result must be identical to the change of 



= 0, 
= 0. 



(4.7) 
(4.8) 
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< 9mn > (0) on the lhs, i.e. to £ m ,„ (0) + £ n , m (0). This results in two further conditions 
which are identical with (3.14), (3.15). 

It is worthwhile to mention that the covariance conditions (4.10) - (4.15) were fulfilled 
in our original approach which was mentioned in the motivation at the beginning of sect. 
3. On the other hand, (3.14), (3.15) failed. This failure is due to the fact that there 
remained contributions which were not averaged. 

We have found an averaging formula which is covariant. It contains 11 constant 
parameters A^, • ■ ■ , (3 multiplying /(r), 6 multiplying F(r)/r, and 2 multiplying 
G(r)), and has to fulfill 8 independent conditions (3.14), (3.15), (4.10) - (4.15). There 
is still some freedom which one can use. For reasons to become clear in the next section 
we impose three further conditions, 

C [a] =D [a] , iora = f,F/r,G. (4.16) 
This implies a further symmetry of K^ n : 



Kmnkl — Kklmn- (4-17) 

The terms with A, B, E, F obviously respect this symmetry automatically. 
The conditions now fix the parameters uniquely. The result is 





A 


B 


c 


D 


E 


F 


f 


3/16 


-21/16 


21/16 


21/16 


-3/16 


-15/16 


F/r 


-21/16 


27/16 


-63/16 


-63/16 


9/16 


225/16 


G 

















-15/2 



5 Covariant averaging of the Einstein tensor 

Besides the metric, the Einstein tensor is the most important object in general relativity, 
because it enters, together with the energy momentum tensor, directly the field equations. 
It would be highly desirable if one could average the Einstein tensor in exactly the 
same way as the metric tensor, and if the averaged Einstein tensor would be identical 
to the Einstein tensor derived from the averaged metric. The old problem that the 
averaged equations are not identical to the equations with the averaged metric would 
then disappear. 

Let us thus investigate the Einstein tensor. In first order of the perturbation one has 

9(1 — W 4-h i — W ■ —h { ■ —h l i S -Uh-^S 

Indices are raised and lowered with Sij here, so their position is in fact irrelevant. The 
averaging formula (3.4) is now used for an arbitrary point x, the integration variables 
are denoted by a prime, and K™ n (x' — x) depends on the difference x' — x. The distance 
r' now means r' = \x' — x\. It is convenient to introduce the modified expression 

K k ±(x> -x) = ^K k ±(x> - x), (5.2) 
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such that 



J K k ±(x> - x)h kl (x')dr>— = j K k ±(x> - x)h M (x')d 3 x'. (5.3) 

This makes partial integrations easy in the case under consideration. Differential oper- 
ators d/dx' n acting on h k i(x') can be shifted to Kfj(x' — x) by partial integration, and 
finally be replaced by d/dx n acting on K^(x' — x). 
There are now two possibilities of averaging: 

The first possibility is to average the metric in the way described before. Subsequently 
one calculates the Einstein tensor from (5.1), using the averaged metric on the rhs. The 
factor of hki in the integrand then becomes 

jkl j/~ikl _|_ f/-kl i jsikl j/~ikl _ T^ikl j r , fy-kl ij x /r a\ 

1 mn ~ 21 i imn ^^mnii 21 m mi J^ n imi J^i ij Omn "T , mn . K -^) 

The second possibility is to calculate the Einstein tensor G kt within the old metric and 
then average it with our formula in exactly the same way as we averaged the metric 
tensor. Shift the partial derivatives from the metric to K, and rename dummy indices 
where necessary such that h k \ appears in all six terms. The factor of h k i in the integrand 
now becomes 

jkl f>ij ftkl _|_ f>kl i _l(ki I k _f?i j Xkl _i_ f>i kl (r r\ 

°mn 1 ^mnnj u ' ly mnH ly mnH ly mnH ly mni>j u ' ly mni> • \ -°J 

We have to check whether the two expressions are identical. This would be an extremely 
complicated task if attacked by brute force. Fortunately one can simplify the problem a 
little bit, although it stays complicated. An inspection of (5.4), (5.5) shows the following 
property. The terms K™ n ,\ are identical. We next use that the tensor K mnk i is invariant 
under the exchange (mn) <-> (kl). The remaining five terms, if arranged properly (e.g. 
the first in (5.4) and the last in (5.5)), correspond to each other by using this symmetry. 
This implies J mnk i = himn- Therefore the difference l mnk i - J mn ki = Imnki - himn 
is antisymmetric under the exchange (mn) <-> (kl). On the other hand, the tensor 
composition of this expression must have the general form (3.6). The terms ~ A, B, E, F 
are symmetric under the exchange (mn) <-> (kl) and thus cannot appear, therefore the 
difference has to be of the form 

T T ( X ~ X )k(% ~ %)l r ( x — x )m(% — x )n-i 
J-mnkl — Jmnki ~ C(T )[0 mn k l J. [0.0) 

The knowledge of this structure allows a considerable simplification, because one may 
now, e.g. contract k — I in order to extract the function c(r'). One cannot further 
contract m = n because then the rhs of (5.6) vanishes identically. Instead of the four 
indices present originally, one thus has to deal with two indices only and can investigate 
the expression 

Imnk ~ Jmnk = C (, r ')l$mn — 3- ^—^ — ]. (5.7) 

It is now necessary to insert the explicit form of K® n (x' — x), to perform the differenti- 
ations, and to make use of the relations between the functions f(r),F(r)/r,G(r). The 
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elementary but tedious calculation gives the structure (5.7) with c(r) a superposition of 
five terms which are proportional to f"(r), f'(r)/r, f(r)/r 2 , F(r)/r 3 ,G(r)/r 2 . 

Let us, for the moment, keep the parameters free and only make use of the special 
forms of the tensors KM, n and K^ mn in (3.12), (3.13), as well, of course, of the con- 
ditions C[ ] = D[ a ] in (4.16). The result is striking, therefore we show the explicit result 
here: 

C ( r ) = ./"(r) 

+ [2A [f] + 2B [f] - B [F/r] - C[F/ r ]] f(r)/r 

+ [-AB [f] + 7B [F/r] + 6C [F/r] + F [F/r] + F [G] ] f{r)/r 2 

+ [-15B [F/r] - 10C [F/r] - F [F/r] ] F(r)/r 3 

+0-G(r)/r 2 . (5.8) 

A sort of miracle happens. All five coefficients in front of the functions vanish for the 
parameters in (4.18)! This implies that the averaged Einstein tensor is identical to the 
Einstein tensor of the averaged metric, i.e. the averaged field equations are identical to 
the field equations of the averaged metric. 

This property is highly welcome and it is hard to believe that it is accidental. Clearly 
the symmetry relations shared by our averaging formula and by the Einstein tensor 
played a central role in the derivation of this result. This becomes evident if one writes 

2G rnn = T™ n h kl , (5.9) 

with the operator 

-5 mn 5 kl & i d J +5 mn d k d l . (5.10) 

Both expressions, K mnkl as well as T mnkh are symmetric under m <-> n, under k <-> I, and 
under (m, n) <-> (k,l). These symmetries implied the vanishing of the five symmetric 
tensors in the difference of (5.4) and (5.5). The vanishing of the remaining antisymmetric 
structure in (5.6) could be demonstrated explicitly, but at present we are not aware of 
some deeper reason behind this. 

The result for the covariant averaging of the Einstein tensor is certainly not trivial. 
For the Ricci tensor, which does not fulfill the above symmetry properties, the relation 
is not valid. 

6 Static perturbations in Minkowski space 

Our extension to the four dimensional case is rather modest. We assume that one can 
find a system in which the perturbation is approximately static. We also neglect the slow 
time dependence in the Robertson Walker metric, so one may choose coordinates such 
that the unperturbed metric is the Minkowski one, (-1,1,1,1). To keep this situation, 
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only rigid translations, rigid spatial rotations, and infinitesimal transformations which 
keep the time unchanged are allowed. This means that £° = 0, and £ m is independent of 
t. Furthermore we can drop all time derivatives in the metric. Under these restrictions 
the perturbations h 00 and h m0 become gauge invariant. 

We average the perturbation with the following simple ansatz. 

</w>(0) = Jk k ±h kl r 2 dr^, (6.1) 

<h 00 >(0) = J p(r)h 00 r 2 dr^, (6.2) 

< h m0 > (0) = J q( r )h m0 r 2 dr^. (6.3) 

Of course the functions p(r) and g(r) have to fulfill the normalization conditions 

/ p(r)r 2 dr = / q(r)r 2 dr = 1. (6.4) 
Jo Jo 

The equations (6.1) - (6.3) are written in such a way that the volume element d 3 x/4n 
appears in the integrals. Hopefully, the fact that we sometimes include the factor r 2 of 
the volume element into the function (as in f(r)), and sometimes don't (as in p(r), q(r)), 
does not produce too much confusion. The choice is motivated by the way how we have 
to perform the various partial integrations. 

Equations (6.1) - (6.3) are a simple generalization of our previous formula. For static 
perturbations they are covariant in the sense of (1.1) with respect to static transforma- 
tions. One could use a more general ansatz, where a term with is inserted into the rhs 
of the averaging formula (6.1), and a gauge invariant combination of h\. and (x k x l /r 2 )hki 
into (6.2). We found that this does not help to solve the problem which will arise at the 
end, therefore we keep things simple and work with (6.1) - (6.3). 

Let us now consider the (lowest order) Einstein tensor which reads 



^^mn — z,u mn r "Oimn l <'0H u mn> 

2Goo = /',,•'' (6.6) 

2G m o = h m o,l —hio, l m ■ (6.7) 

Here G$ n is the spatial part of the Einstein tensor in (5.1). Under the assumptions above, 
the additional terms in G mn , as well as Goo an d G m o are invariant under infinitesimal 
static transformations. 

We now investigate whether the Einstein tensor of the averaged metric can be iden- 
tical to the averaged Einstein tensor. For G m0 this is trivial, we start with G mn . For the 
part G$ n we know from the previous section that the result is independent of the order 
of averaging. We can restrict to the additional contributions in (6.5). 

If we average the perturbation Hq according to (6.2) and introduce into (6.5) we 
obtain the integrand 
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\p,mn -P,l S mn ]h° = [(p" - t.p!^L - (p» + V -)5 mn }hl (6.8) 

ry ry £, ry 

If, alternatively, we first calculate 2G mn in (6.5) with the old metric and then average it 
in the same way as (6.1), i.e. replace hki by Gki there on the rhs, and shift the partial 
derivatives from the perturbation to the multiplying functions, we obtain the integrand 

[K k ±, H -K k mnk ,\K (6.9) 

This gives the condition 

KLM ~K k mnk !l = (f - - (P" + jVmn. (6.10) 

We next apply the same procedure to G o- If we average the perturbations h\ and 
according to (6.1) and introduce into (6.6) we obtain the integrand 

[kPi-KgWhu. (6.11) 

If, alternatively, we first calculate 2Goo in (6.6) with the old metric and then average it 
in the same way as (6.2), i.e. replace /iqo by Goo there, and shift the partial derivatives 
from the perturbation to the multiplying functions, we obtain the integrand 

lb] I 

{p>5 kl -p kl }h kl = -[(p" - ?-)^f - (p» + P -)5 kl \h kl . (6.12) 

■/ ry ry Zj ry 

The condition Kf\) -K k J, ij = -[(p" - ^)^- - (p" + ^)5 kl ] which arises from (6.11), 
(6.12) is a consequence of (6.10) if one renames the dummy indices k, I there, subsequently 
replaces to, n by k, I, and uses the symmetry K mn ui = K k i mn . We thus only need to 
consider (6.10) in the following. 

With (3.5), (3.6) and the parameters in (4.18), a lengthy but elementary calculation 
gives 

K k ±,M -K k mnk n = KArf-^ + K 2 (r)5 mn , (6.13) 

with 



Kl(r) = Am + ^m + ?.m. 30 m t (6 , 4) 

4 t 2 r 2 r r 



W = im-^-^ + K<M. (0,5) 



If one equates the factors of x m x n /r 2 and 5 mn in (6.10) one obtains the two conditions 



K^r) = p"(r) -p'(r)/r, (6.16) 
K 2 (r) = -p"{r)-p'(r)/r. (6.17) 



The sum of these equations gives 
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p , (r) = ~(K 1 (r) + K 2 (r)). (6.18) 

If one calculates p"{r) from this, and inserts again into (6.16), (6.17) one obtains the 
integrability condition 

2K, (r\ 

K' 1 (r) + K' 2 (r) + ^-t = 0. (6.19) 

r 

This condition consists, in fact, of six conditions, i.e. all the constants in front of 
f"'(r)/r 2 , f"(r)/r 3 , f(r)/r 4 , f(r)/r 5 , F(r)/r 6 ,G(r)/r 5 have to vanish. Again a sort of 
miracle occurs. All these constants vanish for the parameters in (4.18). 

Up to now everything worked perfectly, but now we run in trouble with the normal- 
ization condition (6.4). From (6.18) one gets 



I f°° rT 3 1 co 1 f'°° 

1=1 r 2 p(r)dr = —p(r) — - / r 3 p'(r)dr 
Jo ' L 3 J o 3 Jo 

1 r°° 



K 2 {r))dr = -. 



(6.20) 



Unfortunately the normalization becomes 3/2 and not 1. The boundary values [r^r)]^ 
cannot help, because these have to vanish in order that j °° r 2 p(r)dr converges. 

Thus we "almost" succeeded to extend the applicability of the covariant averaging 
formulae (6.1) - (6.3) to the Einstein tensor, but we failed at the end. There is a clash 
between the condition from covariance and the normalization. A slightly more general 
approach can probably solve this problem. 



7 Outlook and conclusions 

The covariant averaging procedure presented in this paper is complicated. This was to be 
anticipated, one could not expect to obtain a simple solution to a complicated problem. 

Compared to most approaches in the literature our formula shows the following fea- 
tures. First one may be surprised that there is no factor ^fg in the integrand. We have no 
comment on this, it simply is not present. An important point is, that it is not sufficient 
to work with bivectors which mix the indices. It is necessary to have bitensors in (3.6). 
Our formula is more general than prescriptions which only average over a certain volume 
because we can use an arbitrary normalized smearing function /(r). But it is essential 
that, besides the function /(r), also the integrals F(r) and G(r) appear. Two simple 
examples suggest themselves. We recall that f(r) should behave ~ r 2 for small r. 
Averaging over a sphere of radius r : 

f(r) = — e(ro-r), F(r) = (-1 + -)6(r -r), G(r) = (- — + _)6(r -r). (7.1) 

r r zr ^ r 

Averaging with an exponentially decreasing function: 

f(r) = ^e-/-, F(r) = -(1 + L + ^K^, G(r) = -(^ + ^K^ . (7.2) 
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One has to face the fact that there is no chance to find a much simpler covariant averaging 
formula in three dimensions than the one presented here. If one works in first order of the 
perturbed metric, there must be a linear connection between the original perturbation 
hki(x) and the averaged < h mn > (0), represented by a tensor (tensor in the sense of 
linear algebra) K^ n which is symmetric with respect tom«ri and to k I. Such a 
connection has also been discussed by Boersma [8 J . The only objects which are available 
for the construction of K® n are the vector x and Kronecker deltas. Therefore one ends 
up with the six tensors defined in (3.6). The functions in front of the tensors must 
depend on r only. They are related by covariance. A special solution has been presented 
in this paper. There are also more general solutions not mentioned here, e.g. averaging 
formulae which, besides f(r),F(r)/r,G(r), also contain the derivative r/'(r). These 
may allow to find a prescription which is also suited for a covariant averaging of the 
Einstein tensor in four dimensions (for static perturbations). The most general form of 
a covariant averaging formula is under investigation. 

A feature of our averaging formula, which might be considered as unpleasant, is the 
appearance of the integrals F(r)/r and G(r), with F(Q) = —1 and G(0) finite. Although 
the term with F(r)/r does not cause a singularity at the origin because the angular 
integration vanishes there, it somehow hampers the smoothing procedure. Before trying 
to find a "better" solution one should, however, take notice of the following fact. It is 
inevitable that functions appear which are not too smooth at the origin. This should 
be clear from our proof of covariance in sect. 4. An infinitesimal transformation within 
the averaging formula must result in the corresponding transformation of the averaged 
metric at the origin. Technically this can only arise through boundary terms at zero 
which originate from partial integrations. Something substantial must be present near 
zero in order to produce these boundary terms. 

An obvious task to be done is to investigate iterations of our averaging formula and to 
check whether the iteration procedure converges. But it is clear that one will not always 
obtain a smooth metric in the limit. The reason is again covariance. We are still free to 
perform gauge transformations, and by an unfavorable choice of gauge the "smoothed" 
metric can look wavy and irregular. All one can expect is that the final metric becomes 
equivalent to a smooth metric. 

One could also proceed to extend the approach to second order in the perturbation. 
We recall that (3.4) also holds if g mn and gui are replaced by h mn and hki- In second 
order an additional contribution which is quadratic in the perturbation will probably be 
needed in the integral. Though certainly tedious, such an extension appears feasible. 

A particularly pleasant property of our formula is the fact that it likewise yields a 
covariant averaging of the Einstein tensor and thus of the energy momentum tensor in 
three dimensions, an extra bonus which, quite surprisingly, came out from the suggested 
averaging formula. We hope to resolve the minor problems found in the four dimensional 
case by a slightly more general prescription. 

Acknowledgement: I thank Juliane Behrend for valuable discussions and for her 
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